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SUMMARY 


We  describe  the  behaviour  of  the  electrons  in  a  gas  discharge  by  using 
the  Boltzmann  equation  and  Fermi  golden  rule  ideas.  The  resulting  equations 
are  applied  to  Sulphur  Hexafluoride  and  electron  energy  distribution  functions 
are  calculated.  The  distributions  are  shown  to  be  in  close  agreement  with 
those  of  a  Monte-Carlo  simulation.  The  ideas  developed  here  should  be  easily 
extended  to  solve  anisotropic  scattering  problems  and  may  also  make  use  of 
the  pseudopotentials  which  are  now  becoming  available  in  the  literature. 
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1  INTRODUCTION 


In  Gas  Discharge  Lasers  the  molecular  energy  levels  are  excited  by  hot  electrons  which 
are  accelerated  to  high  speeds  by  the  applied  electric  field.  In  order  to  appreciate  and 
improve  the  efficiency  of  these  systems  it  is  necessary  to  understand  the  details  of  this 
process.  The  statistical  nature  of  the  scattering  events  means  that  the  electrons  in  the 
electron  swarm  have  a  distribution  of  energies,  in  this  report  we  concern  ourselves  with 
the  evaluation  of  this  distribution. 

There  are  two  established  methods  for  tackling  this  problem.  The  first*,  which  we  will 
consider  in  detail,  introduces  a  distribution  function  at  line  one.  The  distribution,  which 
describes  the  probability  of  finding  an  electron  in  phase  space,  obeys  the  Boltzmann 
transport  equation,  which  may  be  solved  approximately  in  a  simple  regime. 

The  second  method2  treats  the  electrons  individually  and  keeps  track  of  their  phase  space 
trajectories.  This  method  is  exact  and  extremely  useful  for  understanding  spatial  and 
magnetic  field  effects  in  gas  discharges,  however  such  calculations  use  large  amounts  of 
computer  time  and  usually  generate  much  more  information  than  is  needed  for  the  type  of 
problem  which  concerns  experimental  workers  at  the  present  time.  Nevertheless  the 
Monte-Carlo  is  a  valuable  tool  in  so  far  as  it  is  exact  and  serves  as  a  useful  touchstone 
for  Boltzmann  equation  results.  For  the  present  we  concern  ourselves  with  the  Boltzmann 
equation  and  uniform,  steady-state  electron  distributions.  The  Monte-Carlo  method  will  be 
applied  to  dynamic  non-uniform  situations  in  a  later  document. 

In  Section  2  we  introduce  the  electron  distribution  function  and  show  that  it  obeys  the 

Boltzmann  equation.  We  decompose  the  distribution  into  s  and  p-like  behaviour  and  use 
this  to  give  two  coupled  equations. 

In  Section  3  we  evaluate  the  s-type  collision  term  derived  in  Section  2  for  recoil 

scattering  and  inelastic  scattering.  Our  method  is  different  from  the  usual  technique 

employed  in  gas  systems,  we  use  Fermi-Golden  rule  ideas  to  relate  the  transition  rate  to 
the  scattering  cross-section.  Although  our  results  are  the  same  as  those  usually  presented 
we  feel  that  in  view  of  pseudo-potential  data  now  becoming  available  for  gas  atoms'*,  the 
method  is  the  natural  one  to  use. 

In  Section  4  the  p-type  collision  term  is  considered  and  a  momentum  relaxation  time  is 

introduced.  The  equations  derived  in  Section  2  are  then  decoupled.  We  indicate  how  the 
equation  dealing  with  the  s-type  part  of  the  distribution  is  modified  for  a  gas  mixture. 

This  equation  is  the  working  equation  for  the  gas  discharge  system. 

In  Section  5  the  working  equation  derived  in  Section  4  is  solved  under  special  conditions 
to  give  the  Maxwell-Boltzmai.n  and  Druyvestyn  distributions.  In  Section  6  we  describe  two 
general  numerical  methods  for  solving  the  Boltzmann  equation  when  both  recoil  scattering 

and  inelastic  scattering  are  important.  In  Section  7  we  reconsider  the  relaxation  time 
approximation  and  introduce  a  modification  to  the  equations  presented  in  Section  6.  This 
method  is  then  used  to  calculate  the  electron  distribution  function  in  SFg  and  this  is 

briefly  compared  to  the  exact  results  calculated  by  Monte-Carlo  simulations.  For 

completeness  we  consider  distributions  with  the  same  average  energy  and  we  present  a 
graph  showing  the  difference  between  Maxwell-Boltzmann,  Druyvestyn  and  the  SF5 
distribution.  In  Section  8  we  present  a  few  concluding  remarks. 
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2  THE  BOLTZMANN  EQUATION 


We  start  by  introducing  a  distribution  function  f(k,  r,  t)  which  is  defined  as  the 
probability  of  finding  an  electron  with  wavevector  k  at  position  r  at  time  t  with  spin  up. 
The  growth  of  f(k,  r,  t)  is  determined  by  the  current  of  electrons  into  the 

six-dimensional  phase  space  element.  The  electrons  move  in  r  space  by  virtue  of  their 
velocity,  whilst  motion  in  momentum  space  is  due  to  the  force  acting  on  them.  In  the 
absence  of  scattering  we  equate  the  growth  of  f(k,  r,  t)  to  the  divergence  of  this  six 
dimensional  current. 


df  (k.r.t)  -  -  V  .  (fv)  -  V  k.  (fF)  -  -  v  .  Vf  -  F.VRf  , 
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In  the  absence  of  applied  fields  and  concentration  gradients  —  (k,r,t)  is 
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determined  solely  by  collisions 


9f 

ST 


(k,r,t) 


8f 

3F 


(k.r.t) 


col  1 


(3) 


The  Boltzmann  equation  combines  these  two  processes, 
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We  shall  examine  the  simplest  case  when  B  =  0  and  spatial  gradients  are  unimportant. 
We  shall  also  be  concerned  with  the  steady  state.  Then  from  equations  4  and  2  we  have 


—  E  .  V  f(k) 


col  1 


(5) 


where  f  (k)  is  now  a  function  of  k  alone.  Because  the  momentum-transfer  cross-section 
is  large  and  the  scattering  event  is  velocity  randomising  f(k)  will  tend  to  be  isotropic, 
however  we  know  that  the  electron  swarm  must  drift  under  the  applied  electric  field  so 
there  must  also  be  an  anisotropic  contribution  to  f  which  weights  the  distribution  in  the 
direction  of  the  electrostatic  force.  In  general  f(k)  can  be  expanded  in  terms  of  the 
spherical  harmonics 
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y  f  (k)  ?n  (cos  *), 
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where  \p  is  the  angle  between  the  electric  field  and  the  k  vector  direction.  Y/e 
only  the  first  two  terms  are  significant  (the  Lorentz  approximation).  Hence, 


f  (k)  -  f  (k)+  cos  $  f  (k)  . 
o  1 

on  inserting  7  into  5  we  see 
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where  the  axes  have  been  chosen  so  that  k2  is  in  the  direction  of  E,  (see  Fig  1). 
Coverting  to  polars  we  write, 
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Multiplying  equation  10  by  sin  ^  and  integrating  from  0 


to  7r  gives 
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Multiplying  equation  10  by  cos  ^  and  integrating  from  0  to  x  gives 
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We  see  that  we  have  separated  equation  10  into  two  equations  by  considering  s  and 
p-type  behaviours.  The  collision  terms  still  need  to  be  considered  in  some  detail  before 
the  coupled  equations  11  and  12  can  be  solved. 


There  are  two  contributions  to 
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col  1 


the  growth  of  f(k)  is  due  to  transitions  from  all  other  k’  states  into  k,  minus  a  term  due 
to  electrons  leaving  k. 


So, 
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KX 
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where  P  (k\  k)  is  the  probability  per  unit  time  of  a  transition  from  a  full  state  at  k'  to 
an  empty  state  at  k.  The  V/8z^  factor  is  the  density  of  states  in  k  space. 

Substituting  equation  7  into  13  we  obtain 
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We  assume  that  P(k,k')  depends  on  k,k*  and  the  scattering  angle  (7)  and  not  on  the 
absolute  orientation  of  the  vectors  in  space. 

It  is,  therefore  convenient  to  evaluate  the  d^k'  integrals  in  a  frame  which  has  the  kz 
directed  along  the  initial  k  direction.  In  this  case 


cosi£'  —  cos cosy  +  sin0  siny  cos/3  ,  (15) 

where  j3  is  the  azimuthal  angle,  (for  details  see  Holtstein^) 

Substituting  (15)  into  (14)  and  carrying  out  the  integrals  to  separate  the  two  contributions 
we  obtain 
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These  are  still  formal  equations,  in  the  next  section  we  look  at  the  influence  of  recoil 
scattering  and  inelastic  scattering  on  equation  17  and  in  the  following  section  we  examine 
equation  16. 

3  COLLISION  PROCESSES 

So  far  we  have  said  nothing  about  the  details  of  the  microscopic  events  which  remove 
energy  from  the  electron  distribution.  These  fall  into  two  main  categories,  recoil  and 
inelastic  scattering.  Under  normal  circumstances  recoil  scattering  has  little  direct  effect  on 
energy  relaxation  as  only  a  small  amount  of  energy  is  removed  in  each  scattering  event, 
however  it  can  be  important  at  low  electric  fields  when  the  second  process,  molecular 
excitation  is  unimportant  because  the  electrons  have  not  yet  got  enough  energy. 

In  this  section  we  evaluate  the  s-type  part-  of  the  collision  term  by  considering  a 
Fermi-Golden  Rule  approach.  Smith  and  Thompson*  used  hand  waving  arguments  in  their 
derivation  and  rigorous  derivations  seem  in  the  main  to  be  absent  from  the  literature 
although  Holstein’s^  method  is  correct.  Our  method  relates  the  transition  rate  to  the 
potential  responsible  for  the  scattering  and  in  view  of  the  pseudo-potentials  now  available** 
in  the  literature  for  calculating  scattering  rates  by  just  such  a  route  it  seems  like  a  worthy 
one  to  pursue.  We  start  by  quoting  Fermi’s  Golden  rule:- 
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Given  an  electron  in  state  k,  the  transition  rate  to  an  empty  state  k'  is  given  by  Fermi's 
golden  rule. 


2z 


P(kfk')  -  -  <k  eV  (r)  k‘>  o  (e  (k)  -  e  (k ')  -  A  (k.k')), 


(18) 


where  A(k,k')  is  the  energy  tost  by  the  electron  in  the  transition  and  V(r)  is  the 
perturbing  potential  due  to  the  atom. 

We  now  examine  the  two  types  of  scattering  process. 


3.1  Recoil  Scattering 

We  assume  for  the  moment  that  all  the  atoms  are  at  rest  before  a  collision  happens  ie 
temperature  =  0.  Any  energy  dumped  into  the  molecular  system  is  also  assumed  to  be 
removed  instantaneously.  So  energy  from  the  molecular  system  may  not  feed  into  the 
electron  system.  The  scattering  rate  is  usually  written  as 


r  =  nQv, 


(19) 


where  Q  is  the  cross-section,  v  is  the  speed  of  the  electron  and  n  is  the  number  of 
scatterers  per  unit  volume.  If  the  nature  of  V(r)  is  known  we  may  also  evaluate  T(k) 
through  the  relation 


r  (k)  - 


P  (k,k‘ ) 


(20) 


and  equation  18. 


Now  for  electrons  with  an  energy  less  than  30  eV  we  shall  assume  that  the  scattering 
centre  is  much  smaller  than  the  wavelength  of  the  electron.  (Typically  only  ~  10  times 
smaller).  This  allows  us  to  replace  the  potential  V  (r)  with 

V(r)  -  S06(r-r0) ,  (21) 

where  S0  is  the  scattering  strength. 

If  we  treat  the  collision  as  elastic  then  we  may  easily  evalulate  the  scattering  rate  T(k) 
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where  m  is  the  mass  of  the  elctron  and  we  have  used  equations  18,  20  and  21. 
So  we  see 


(23) 


d3r  (24) 


If  V(r)  is  finite  in  extent  then  it  can  be  shown  by  partial  wave  analysis  that  equations  23 
and  24  still  hold  approximately  if  the  spread  of  V(r)  is  much  smaller  than  the  electron 
wavelength^.  Experimentally  however  Q  is  dependent  on  energy.  This  dependence  is  due 
to  the  failure  of  our  sharp  potential  assumption.  Indeed  the  Ramsauer-Townsend  minima 
in  experimentally  measured  Q(e)  values  is  a  manifestion  of  resonant  tunnelling  through  the 
potential  barrier  due  to  the  size  of  the  atonA  To  include  the  energy  dependence  in  a 
simple  way  and  still  retain  the  isotropic  scattering  ideas  we  write  the  transition  rate  as 


where 


V  (r) 


2z 

P ( k ’ ,k)  -  —  A(f’)  b  (('-(  -A)  ,  (25) 

* 

where  A  is  the  energy  lost  in  a  transition  from  k'  to  k  and  A(f')  is  a  function  of  the 
initial  energy  of  the  electron.  Then 


A(0 


Q  ( ( )  t 
V  m2 


(26) 


We  have  retained  the  isotropic  scattering  approximation  whilst  incorporating  the  energy 
dependence  of  Q(().  1°  reality  anisotropic  scattering  and  the  energy  dependence  of  Q(r) 
are  inextricably  linked  to  the  failure  of  the  notion  of  delta  function  scattering  centres.  As 
we  are  mainly  interested  in  the  energy  dependence  of  the  electron  distribution  function  it 
is  hoped  that  such  an  approximation  will  not  lead  to  large  errors  in  calculated  excitation 
rates. 

In  view  of  our  preliminary  remarks  we  note  that  recoil  collisions  can  only  remove  energy 
from  the  electron  system.  Hence  the  first  term  in  equation  17  which  is  the  scattering  in 
term  (scattering  into  k)  must  have  k  <  k\  similarly  the  second  term  has  k'  <  k  as 
electrons  can  only  be  removed  from  k  by  a  decrease  in  energy. 

For  recoil  scattering  the  initial  and  final  energies  are  related  through  the  equation  1 
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where  e*  is  the  energy  of  the  initial  state  k\  (  that  of  the  final  state  k,  7  is  the 
scattering  angle  and  M  is  the  mass  of  the  recoiling  atom. 

Defining 
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we  may  write 


P(k\k)  -  —  A ( € '  )  5  (g  (f\f))  (29) 

We  may  now  re-write  equation  17  solely  in  terms  of  energy 
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we  note  that  g  (e',e)  =  0  is  approximately  satisfied  by  e  =  c'. 
We  use  the  property  of  the  delta  function 
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where  the  xn  are  the  set  of  numbers  which  satisfied  g(x)  =  0. 

Now, 

I  g'  («',  O  I  ~  1  -  -  (1  -  cos  7)  (32) 

1  1  M 


and 


g'  (f ,  € ' )  I  ~  1  +  -  (1  -  COS  7) 

1  M 


(33) 


8 


Hence  we  obtain, 
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now  £+A~$m(v  +  av)2 

where  v  is  the  speed  of  an  electron  with  energy  e  and 


a  «  ~  (1  -  cos  7) .  (35) 

Writing  everything  in  terms  of  speed, 
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where  we  have  integrated  over  0. 
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Expanding,  integrating  and  using  equation  26  we  find 
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converting  to  an  energy  derivative, 
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This  result  aloig  with  the  equation  derived  in  the  next  section  will  allow  us  to  complete 
equation  1 1 . 


3.2  Inelastic  Scattering 

As  well  as  being  able  to  reduce  its  energy  by  increasing  the  kinetic  energy  of  the 
surrounding  gas  molecules  an  electron  may  excite  the  gas  molecule  by  changing  either  its 
vibrational  state  or  by  forcing  an  electron  transition.  This  is  usually  the  dominant  process 
for  energy  relexation  in  a  typical  discharge  as  these  transitions  relax  large  amounts  of 
energy  (typically  ~10  eV).  In  this  section  we  consider  how  this  scattering  mechanism 
influences  the  s-type  part  of  the  collision  term. 

We  shall  assume  that  the  transition  rate  may  again  be  written  in  the  form 


P  (k,k' )  -  —  A(()  6  («-«*  -  4), 
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(39) 


where  A  is  now  a  constant  corresponding  to  an  excitation  energy  in  the  gas  molecule  and 
P(k.k')  takes  an  electron  to  a  reduced  energy.  Evaluation  of  the  scattering  rate  T(k)  gives 
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equating  T  with  NQv  we  find 
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We  now  examine  the  8f  /8t  |  ,,  due  to  inelastic  collisions 
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substituting  equation  41  into  42  gives 
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f o ( k ' )  A  (e+A)  J  e+A  -  fQ(k)  A(e)  J  e-A 


2m  ~  (43) 

*4 


2  2 
jj  k> 

where  k'  is  now  the  solution  of  -  -  e  +  A 


2m 


substituting  (41)  into  (43)  gives 


8f 

0 

3T 
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f  (e+A)  NQ  (e+A)  (e+A)  -  f  (e)  NQ(e)  e 
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inel 

V  v 

(45) 


mf 


This  result  completes  equation  11,  however  before  the  coupled  equations  11  and  12  can 

3f, 


be  solved,  we  need  to  calculate 


1 


8t 


col  1 
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4  THE  RELAXATION  TIME  APPROXIMATION 


In  the  presence  of  an  applied  electric  field  the  distribution  function  f(k)  contains  an 
anisotropic  contribution  fj(k)cos0  which  would  decay  if  the  field  was  turned  off.  We  model 
this  decay  with  the  ansatz. 


8fl 

—  (k) 

3t 


col  1 


fj  (k) 
t  (k) 


(46) 


where  r(k)  is  a  momentum  relaxation  time  which  is  a  function  of  k  alone. 

If  we  assume  that  the  main  contribution  to  the  relaxation  rate  is  due  to  quasi-elastic 
scattering  through  the  recoil  scattering  process  then 


P(k.k’)  ~  P(k\k)  and  f  (k’>  *  f  (k) 
Substituting  46  into  16  we  find, 


f1  (k) 
r  (k) 


(cosy  -  1)  fj(k)  P  (k*  ,k)  d3k’ , 


(47) 


which  gives  the  momentum  relaxation  rate, 


1 

t  (k) 


(1-cos  7)  P  (k* ,  k)  d3k' 


(48) 


For  isotropic  scattering  the  momentum  relaxation  rate  equation  48  and  the  scattering  rate 
equation  20  are  identical.  In  this  case 


i  -  r  -  NQv  (49) 

7 


In  general  this  expression  will  only  be  a  valid  representation  of  the  momentum  decay 
process  when  the  scattering  rate  due  to  recoil  scattering  is  much  larger  than  that  due  to 
inelastic  cross-sections.  Many  authors  have  considered  the  breakdown  of  Boltzmann 
calculations  in  a  gas  to  be  due  to  the  truncation  of  the  series,  however  the  validity  of  the 
relaxation  time  approximation  has  not  been  investigated.  We  shall  assume  that  we  are 
working  at  an  electric  field  where  equation  49  is  valid,  then  equation  12  becomes 

ar 

—  E  — -  -  -  f.  NQv  (50) 

1)  dk 
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Similarly  for  equation  11 
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(51) 


or  in  terms  of  energy, 
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Hence  from  equations  38  and  45 
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4NQ  c2  r 
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3c 
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(53) 


Making  the  change  of  variables  c  “  eu,  we  obtain 


E_ 
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3r 
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2m 
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NQ  u2  f 

NQ 

3u 

M 

3u 

0 
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f  (u  +  A)  NQ  (u  +  A)  (u  +  A)  -  f  (u)  NQ  (u)  u 

O  0 


0  . 


(54) 


This  is  the  working  equation  for  the  single  component  gas  when  the  temperature  is  small 
and  the  relaxation  time  approximation  is  valid 

We  are  usually  concerned  with  a  gas  mixture  at  finite  temperatures,  then  the  equation  is 
modified  to  give^ 


13 


where  k  labels  the  gas  species  and  j  corresponds  to  the  j-th  excitational  state,  whilst 
Qjf(u)  is  the  corresponding  cross-section.  Njj  is  the  total  number  of  k-type  molecules 
whilst  Ng  is  the  number  of  k-type  molecules  in  the  ground  state  and  Nj[  is  the  number 
of  k-type  molecules  in  the  j-th  excited  state. 

Before  we  proceed  to  a  full  numerical  solution  of  equation  55  we  examine  an  analytical 
limit. 


5  THE  MAXWELL-BOLTZMANN  AND  DRUYVESTYN  DISTRIBUTIONS 

There  are  two  well  established  distributions  which  may  be  used  as  best  first  guesses  to  the 
electron  distribution  when  some  quantity  such  as  the  average  energy  is  known  and  others 
need  to  be  estimated.  The  Maxwell-Boltzmann  distribution  is  the  well  known  thermal 
equilibrium  distribution.  The  Druyvestyn  distribution  is  less  well  known  and  is  relevant  to 
the  high  field  regime,  when  thermal  effects  are  less  important.  To  derive  both  and 
consider  the  transition  from  one  to  the  other  we  consider  a  single  component  gas  when 
inelastic  scattering  is  unimportant  and  the  momentum  transfer  cross-section  is  independent 
of  energy,  then  the  Boltzmann  equation  reads. 
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integrating  twice  we  find 


(56) 
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where 
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6mQm 


V 

and  0  -  —  • 


(58) 


and  we  have  dropped  the  k-label  as  it  is  now  superfluous.  Under  weak  electric  field 
conditions,  when  alpha  is  small,  the  shape  of  the  distribution  is  dominated  by  the 
temperature  of  the  gas  molecules,  giving  a  Maxwell-Boltzmann  distribution, 


fQ  -  exp 


(59) 


In  the  field  dominated  regime  we  may  consider  beta  to  be  small,  giving  the  Druyvestyn 
distribution, 


fo  -  exp 


(60) 


The  characteristic  energy  of  the  electron  gas  (eV)  is  given  by  the  ratio  of  the  diffusion 
coeffcient  to  the  mobility. 
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(62) 


On  evaluating  the  integrals,  we  find  that  in  the  temperature  dominated  regime 


€ 


k 


V 


e 


(63) 


which  is  a  restatement  of  the  Einstein  relation  for  a  non-degenerate  electron  gas. 


In  the  field  dominated  regime  we  find 


(64) 


In  general  the  momentum  transfer  cross-section  is  a  function  of  energy  and  inelastic 
scattering  is  important.  In  the  next  two  sections  we  describe  how  general  solutions  may  be 
found  for  a  gas  mixture. 


6.1  CALCULATING  THE  ELECTRON  DISTRIBUTION  BY  GAUSS-SEIDEL  ITERATION 

We  consider  equation  55  and  make  the  transformations  n(u)  =  tju£  fo(u),  where  rj  is  the 
normalisation.  After  some  manipulation  we  find 
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j  k 


where  is  the  energy  needed  to  take  a  k-type  species  into  its  j-th  excited  state, 
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qR  (<  +  (  >  .•  (<)  -  Qk  (<)  ,  (O. 


i*(e)  is  the  velocity  of  the  electron  and  is  the  proportion  of  k-type  gas  molecules  in 
excited  state  j. 


We  discretise  by  using  the  approximate  derivative  expressions 


9n  np+l  np~] 
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where  At  is  chosen  as  the  fundamental  energy  step.  After  some  manipulation  we  find 
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(69) 


The  ak's  and  bk's  have  physical  interpretations  as  the  rates  at  which  electrons  are 
upscattered  and  downscattered  respectively.  It  follows  from  the  definition  of  n(f')  and  the 
finite  nature  of  fo(0)  that  bj  =  0  and  a0  a  0.  Hence  if  the  discrete  Boltzmann  equation 
is  written  down  in  matrix  form  then  to  ensure  that  electrons  are  not  scattered  to  negative 
energies  then  certain  of  the  elements  are  reduced  in  form. 

In  the  absence  of  inelastic  scattering  the  matrix  is  tridiagonal  and  a  simple  analytic 
recursion  relation  exists  for  np.  With  the  boundary  condition  we  find 


^2n2 


alnl 


and 
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the  equation's  are  linear  in  np  and  the  normalisation  is  still  yet  to  be  fixed  we  can 
choose  nj,  then  the  recursion  relation  may  be  used  to  solve  for  np  for  elastic  scattering 
given  any  form  for  Q^(e'). 
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In  general  when  excitation  events  are  important  solutions  of  equation  67  may  be  found  by 
Gauss-Seidel  iteration. 

6.2  DIRECT  SOLUTION  BY  ITERATION 

The  method  quoted  above  can  consume  large  amounts  of  computer  time  before 
convergence  is  achieved.  In  this  section  we  consider  a  method  which  is  quick  when 
superelastic  collisions  are  unimportant.  We  consider  a  direct  solution  of  equation  55. 

Equation  55  may  be  approximated  by  centred  finite  differences^,  with  an  energy  step  of 
du,  to  give. 
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and  6 is  the  nearest  integer  to  uj^/du.  We  have  also  detailed  the  following  quantities 
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In  the  absence  of  superelastic  collisions  =  0,  then  equation  71  may  easily  be  solved. 
A  guess  is  made  for  the  average  energy  of  the  electrons  in  the  distribution.  Typically  this 
will  be  less  than  the  threshold  energy  of  the  first  inelastic  cross-section.  The  maximum 

energy  considered  in  this  calculation  is  set  at  two  or  three  times  this  value,  fj  .  and 

fjfllax+1  are  Set  t0  ^'nite  but  sma"  va'ues  ancl  tJie  difference  is  chosen  sucli*  that 

fImax  flmax+1  as  we  exPect  to  be  in  the  ta‘i  °f  the  distribution.  Values  of 

[fNJN>I  max+1  are  assumed  t0  be  sufficiently  small  so  that  electrons  scattering  down  from 
higher  energies  have  a  tiny  effect  on  the  distribution  function.  Equation  71  is  then 
iterated  down  to  fj. 

This  technique  is  much  quicker  than  the  Gauss-seidel  method  as  the  operation  described 
above  only  has  to  be  carried  out  lmax  times  whereas  the  Gauss-seidel  technique  scales  as 
2 

^max- 

7  APPLICATION  TO  A  REAL  GAS  fSF$) 

The  momentum  relaxation  time  quoted  in  section  4  is  only  approximately  correct,  in 
general  inelastic  scattering  will  affect  it  in  a  none  trivial  way. 

It  has  been  suggested  that  the  momentum  relaxation  rate  should  be  proportional  to  the 
total  cross-section  rather  than  the  momentum  transfer  cross-section.  This  approximation  is 
incorrect  in  detail,  ideally  we  should  find  a  momentum  relaxation  time  which  is  the 
self-consistent  solution  of  equation  16.  In  the  case  of  recoil  scattering  the  quasi-elasticity 
of  the  process  allows  us  to  divide  through  by  the  p-like  part  of  the  distribution  and 
remove  it  from  the  equation.  Inelasticity  on  the  other  hand  allows  no  such  simplification 
and  solutions  of  equation  16  should  ideally  be  found  by  iteration. 
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However  we  have  adopted  the  total  cross-section  method  as  a  best  first  approach  and 
treated  the  ionisation  and  attachment  processes  just  as  inelastic  processes  which  removed 
energy  from  the  electron  bath.  The  only  effect  of  these  modifications  was  to  adjust  the 
first  bracket  in  equation  55  so  that  the  momentum  transfer  cross-section  was  replaced  by 
the  total  cross-section. 

We  focussed  our  attention  on  SFg  and  taking  the  cross-sections  supplied  by  Itoh^  (See 
Fig  2)  we  used  both  numerical  methods  to  produce  electron  distributions.  Both  techniques 
produced  the  same  result  giving  us  confidence  in  our  computer  code. 

The  E/N  value  was  varied  until  the  electron  creation  rate  was  the  same  as  the 
annihilation  rate. 


ie 


Qion  (°  ef(°  df  = 


Qatt  (<)  ff  (e)  df. 


This  gave  us  the  stable  working  field.  A  Monte-Carlo  program^  was  also  run  for  this  E/N 
value  and  the  results  were  found  to  be  in  close  agreement  (See  Figure  3). 

For  completeness  we  have  compared  the  calculated  SF5  distribution  to  Maxwell-Boltzmann 
and  Druyvestyn  distributions  with  the  same  average  energy,  these  results  are  presented  in 
Figure  4.  We  see  that  the  shape  of  the  true  result  is  more  closely  followed  by  the 
Druyvestyn  distribution  which  falls  off  more  quickly  with  energy  than  than  the  slowly 
decaying  Maxwell-Boltzmann.  We  also  note  that  neither  of  the  analytical  distributions 
reproduce  the  tail  of  the  true  distribution  and  the  Boltzmann  calculation  is  necessary  if  we 
are  to  estimate  the  excitation  rates  of  higher  lying  energy  levels  with  any  accuracy. 


CONCLUSIONS 


We  have  used  Fermi  Golden  rule  ideas  to  derive  the  Boltzmann  transport  equation  for  a 
single  component  gas  at  low  temperatures.  A  modified  Boltzmann  equation  has  been  solved 
within  the  relaxation  time  approximation  to  give  good  agreement  with  an  exact 
Monte-Carlo  simulation  for  SF§.  The  Boltzmann  code  is  simple  in  form  and  may  easily 
be  used  to  consider  gas  mixtures  and  to  evaluate  stable  working  fields.  We  note  however 
that  the  relaxation  time  approximation  used  in  this  report  is  inexact,  it  would  be 
interesting  to  calculate  a  self-consistent  relaxation  time  by  an  iterative  procedure. 

In  the  longer  term  Fermi's  Golden  rule  and  the  gas  pseudo-potentials  which  are  now 
available  (and  shown  to  be  accurate1*)  should  make  it  possible  to  calculate  the  recoil 
cross-section  for  simple  gases. 
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Normalised  Distribution  y* f(u) 


whilst  the  curve  results  from  solving  the  Boltzmann  equation. 
( E/n  =367Td) 


Druyvestyn  distribution  function 


as  the  Boltzmann  result  in 


Maxwell-Boftzmann  distribution 


A  Maxweli-Boltzmann  number  density  distribution  with  the 
same  average  energy  as  the  Boltzmann  result  In  FIGURE  3. 
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